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r*"** . Motivated by a recent measurement of proton-proton elastic scattering observables up to 3.0 GeV, 

1 we investigate the description of those data within models of the nucleon-nucleon (NN) interaction 

' valid above the pion production threshold. In addition to including the well known Delta resonance 

s ! , we incorporate two low-lying N* resonances, the iV*(1440) and the AT* (1535), and study their 

j^jfy influence on pp and np observables for projectile laboratory kinetic energies up to 1.5 GeV. 

< 

I. INTRODUCTION 

Elastic nucleon-nucleon (NN) scattering is a process fundamental to our understanding of the nuclear forces. 
This knowledge forms the basis for a broad range of applications from fcw-nuclcon reactions to heavy ion physics. 
Consequently, many experimental and theoretical studies have been devoted to the subject. Since the measurements 
at, e.g., LAMPF [1, 2], the database has doubled, and global phase shift analyses now extend to 3.0 GeV projectile 
kinetic energy for pp scattering and 1.3 GeV for np scattering [3, 4]. Measurements of polarization observables 
(analyzing power, spin correlation coefficients, depolarization parameters and spin transfer triple coefficients) at 
discrete kinetic energies up to 2.9 GeV were obtained by the NNSaturne II Collaboration [5, 6]. At lower kinetic 

■ energies (0.2 - 0.45 GeV), precise analyzing powers and spin correlation coefficients have been provided by the PINTEX 
[ Collaboration [7-9] using the IUCF Cooler Ring and an internal polarized hydrogen target. At higher kinetic energies 

\Q . (0.45 - 2.5 GeV), precise and internally consistent excitation functions, i.e., differential cross sections [10, 11], analyzing 
C*~) ' powers [12, 13] and spin correlation parameters [14, 15], have recently been measured by the EDDA Collaboration 
qq , using the Cooler Synchrotron COSY and internal targets. These new proton-proton elastic scattering data are included 
f^i ■ in the new global phase shift analysis up to 3.0 GeV kinetic energy with its most recent solution SP07 [3]. 

\ The theoretical description of the nuclear force also has quite a long history. During the last decade a lot of effort 

■ has been directed toward deriving the nuclear force in the framework of effective field theory (EFT), especially in 
\ the form of chiral perturbation theory (%PT). This EFT can be regarded as the low-energy limit of the underlying 

■ theory quantum chromodynamics (QCD) and inherits its spontaneously broken chiral symmetry. As first pointed 
out by Weinberg [16, 17], the nucleon-nucleon potential can be expanded in powers of the xPT expansion parameter 

^ P = P'™* , where p ~ is a typical nucleon three-momentum and A x ~ 1 GeV is the chiral scale. Weinberg's idea 
has since been extended to several orders in the chiral expansion with NN potentials of xPT derived to NNLO in 
Refs. [18-20] and to N 3 LO [21, 22]. At N 3 LO the description of NN data below E lah = 200 MeV is comparable in 
quality to that obtained in the so-called "high-precision" NN potential models [23-25] , provided cutoffs in the range 
500-600 MeV are considered [22]. However, the energy range appropriate for the xPT expansion is determined by the 
expansion parameter P being "small" . Thus energy regimes of ~ 0.5 GeV and higher are not accessible by currently 
employed EFT methods. 

In the later 1980's and early 1990's considerable theoretical effort was directed towards understanding the NN 
interaction in the energy region between 0.4 and 1.0 GeV, which is dominated by the Delta isobar. These approaches 
may roughly be divided into two classes [26] . Firstly, there are the models based on a formulation of genuine unitary 
(NNir) three-particle equations, explicitly including the 7rA~ scattering amplitudes P33 [Delta (A) resonance] and Pn 
(Roper resonance) and usually employing separable potentials. The second approach is based on the extension of 
conventional low-energy potentials into the singlc-pion production region [27-30] . Most of these are based on coupled 
channel calculations containing NN, NA and AA channels, and the single pion production is incorporated through 
the decay of the Delta isobar doorway states. The main differences between the two types of models are the choice of 
kinematics, the description of the width of the Delta resonance, and the choice of NN potential being extrapolated. 
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More recent theoretical work on the NN interaction in the energy domain around 0.5 GeV and higher was driven 
by the need of a better understanding of pion production in different reactions, e.g., in pion photoproduction from 
the deuteron in the Delta region [31, 32] and pion production in heavy-ion reactions [33], measured at GSI. The 
exploratory study of Ref. [34] investigates how well a meson exchange model containing the Delta resonance together 
with a Wood-Saxon type optical potential can describe cross sections and polarization observables in the region 
between 1 and 3 GeV. The authors point out that the polarization observables in pp scattering for energies above 
1 GeV arc badly represented by the model employed. 

In this work we want to concentrate on the structure of the nuclear force as revealed in NN scattering between 
1 and 1.5 GeV projectile kinetic energy, the so-called 2nd resonance region. This regime is dominated by low lying, 
well-isolated nucleon resonances. A better understanding of the properties and parameters of those resonances is now 
being pursued at Excited Baryon Analysis Center (EBAC) [35] based on photoproduction data at ELSA (Bonn), 
JLab (Newport News), MAMI (Mainz), GRAAL (Grenoble), and Spring-8 (Hyogo). We want to study the influence 
of the lowest lying irN resonances on the NN data in the region around 1 GeV, which is slightly beyond the region 
dominated by the Delta resonance. We start from a meson exchange model and introduce the lowest order loop 
diagrams containing the A, Roper, and Su resonances into the model. In Section II we give the underlying interaction 
Lagrangians and in Section III we describe the resonance part of the potential and the parameterization of their width. 
In Section IV we describe our data fitting procedure, and in Section V we discuss our results obtained when successively 
introducing the resonances into the potential model. We summarize and conclude in Section VI. 



II. RELATIVISTIC MESON-EXCHANGE MODEL FOR NN SCATTERING ABOVE THE PION 

THRESHOLD 

For projectile energies in the GeV range it is mandatory to employ a relativistic formulation of NN scattering. A 
suitable starting point would be the four-dimensional Bethc-Salpeter equation [36] , which however is very difficult to 
solve. The only successful solution for NN scattering including Delta isobar (A) degrees of freedom was obtained 
in Ref. [28]. Since three-dimensional scattering equations are more amenable for computations, so-called three- 
dimensional reductions of the Bethe-Salpeter equation have been used [29]. In this work we will use the reduction 
proposed by Thompson [37], which was also used in Ref. [34]. For completeness we give the invariant scattering 
amplitude T as introduced by the Thompson reduction: 

r . M 2 1 

T(p',p,W) = V(p',p) + J dk V(p', k) w w + fe _ ^ T(k, p, W). (2.1) 



Here W is the invariant mass of the NN system, W = 2E p , and E p = \/p 2 + M 2 in the cm. system. Redefining the 
scattering amplitude as well as the potential as 

T(p',p,W) = lLf(p',p,W)¥-, (2.2) 

V(p',p) = ^-V(p',p)^, (2.3) 

leads to a scattering equation which has the form of a standard Lippmann-Schwinger equation with relativistic kinetic 
energies 

T(p',p,W) = V(p',p) + J rfk V (p / , k) w + J _ 2E T(k, p, W), (2.4) 

which can be solved with standard numerical techniques. 

The pseudoscalar pion provides the long-range part of the nuclear force and most of the tensor forces. Shorter 
range contributions of the potential are being represented by the exchange of heavier mesons, the vector mesons p 
and oj, representing the 2ir P-wave resonance and a 3-7T resonance. The to provides most of the short range repulsion 
of the nuclear force. The intermediate attraction is provided by the exchange of a light scalar- isoscalar "er" -meson, 
which has not been directly observed. A recent reanalysis [38] of irir data gave evidence for the existence of a very 
broad resonance at threshold with the quantum numbers of the vacuum, a mass of about 440 MeV, and a width of 
about 272 MeV. For a simple model which still can capture basic features of the NN interaction the exchange of it 
(138 MeV), "cr" (~ 500 MeV), p (769 MeV), and uj (783 MeV) is sufficient. The coupling of these mesons to the 
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nucleon is given by the following interaction Lagrangians 

fp 



C pv = -JlLfyf^d^, (2.5) 



TYlps 

C s = -g s ^ip( s \ (2.6) 

£v = - 9v ^M v) ~ ^^"Xa^l" 1 - dvtpW). (2.7) 

Here M is the nucleon mass, and m a , with a = pv,ps, s, v, denotes the masses of mesons with pseudo- vector, pseudo- 
scalar, scalar and vector character. If the mesons carry isospin, the field should be replaced by r • where 
r ? arc the usual Pauli spin matrices. Very often one-boson-exchange (OBE) models also include the r\ (547 MeV) 
and the ao (980 MeV) [39]. However, if the scalar-isovector ao is included, there is no reason to not also consider the 
scalar-isoscalar fo (980 MeV) with similar mass, as well as the rj' (958 MeV). It is clear that mesons with large masses 
will contribute to the short range part of the TV TV interaction. Thus their effect is very small when considering low 
energies, and there they can in principle be omitted. However, since we want to investigate a meson exchange model 
at laboratory energies around 1 GeV and higher, we include for consistency all strangeness-zero mesons with masses 
below 1 GeV. 

In order to introduce contributions from the nucleon resonances into the potential model, we consider the lowest 
order loop diagrams containing a resonance as an intermediate state. The first diagrams to consider are those with 
a nucleon and a single resonance (A, Pn, Sn) in the intermediate state. As a first step we want to consider here 
only low lying S- and P-wave resonances. These resonances can be on-shell at the NN laboratory kinetic energies 
Tub = 632 MeV for the NA, 1136 MeV for the NPn, and 1381 MeV for the NSn diagrams, respectively. These 
energies are given in Table I together with the pion and eta production thresholds. The energy at which two Delta's 
are on shell in an intermediate state is below the threshold for the NSn intermediate state, therefore the diagram 
with A A intermediate states is also included in our calculation. The diagrams included are shown in Fig. 1. 

Since the Delta is an isospin-3/2 state, only isovector mesons can excite it. The only possible exchange mesons are 
then 7T and p. Note that in order to fulfill the weak unitarity bound (a tot > °ci), diagrams with ir-ir, n-p, and p-p 
exchange must all be included. The Pn resonance has the quantum numbers of the nucleon and decays predominantly 
into Ntt (60-70%) [40] and to a lesser extent to Ntttt. We consider here the exchange of a pion and in addition the 
exchange of a "tr", which mimics (correlated) two-pion exchanges. The decay of the Pn into Np has a probability 
of less than 8% and thus this p exchange is not considered. The Sn resonance is the only known resonance that 
couples strongly to the Nij channels. The dominant decay channels are Ntt and N-q [40]. Thus, for this resonance we 
consider only the exchange of tt and rj. Again, in order to fulfill the weak unitarity bound, each of the diagrams in 
Fig. 1 represents in reality four separate diagrams. 

The interaction Lagrangians for the couplings to the Delta are given by 



C NA n = -f^^T^d^ +H.c. (2.8) 



C 



il^^Y'T^ (dr<pU> - d v ^A + Ha, (2.9) 



NAn — —1 

where tp^ is the Rarita-Schwingcr field [41-43] describing the spin-3/2 Delta-isobar and t/j stands for the nucleon field. 
The operator T acts between isospin-1/2 and isospin-3/2 states and H.c. denotes the Hermitian conjugate. 
The Pn resonance has the quantum numbers of the nucleon with the interaction Lagrangians given by 

Jn N* TV 

Cnn* * = ^feM^+H^, (2-10) 

£nn* Pii <j = -9NN^<7 , ipp 11 'tp<Pa+'H-c- (2-11) 

In contrast to the P\\ resonance, the Sn resonance has negative parity, and the interaction Lagrangians are given by 
[44] 

CnNZ tt = -igNN' Tr i'Stt'T tpVir + H.C, (2.12) 

i>ix ^ll 

jC.NNZ n = -igNN* n V'SiiV , ¥>77 + H.c. (2-13) 

* 1 1 a 1 1 

All the meson-baryon vertices described by the above Lagrangian interactions are modified by form factors of the 
dipole-type 

M(p'-p)') = ( A / + v"!> )°°. ("4) 
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where (p' — p) is the momentum transfer between the two interacting baryons. The parameter A Q is the so-called 
cutoff mass and is determined by fitting NN data. We choose n a = \ for all vertices with the exception of the NAp 
vertex, where n a = 2 is applied. 



III. RESONANCE CONTRIBUTIONS TO THE POTENTIAL 

The nucleon resonances enter the NN potential via the loop diagrams given in Fig. 1 . Schematically, the contribution 
to the potential from these "iterative" diagrams can be written as 

VnN,R!R 2 77} p 7} V Rl R 2 ,NN- (3.1) 



Here W is the invariant mass of the NN system, and E^, Eji 2 the energies E Ri = ym 2 ^. + k 2 of the intermediate 
states, which can be either a nucleon, a A(1232), a iV*(1440) or a 7V*(1535). The transition potentials V£ N NA and 
Vnn aa can be found in Ref. [45] and NA and AA in Ref. [46] . The transition potentials for the two N* 
resonances can be calculated in a straightforward fashion, since they carry isospin 1/2. Since the resonances have a 
finite life time, their mass is modified by the width of the resonance according to 

/i ri = m Ti — i—^-- (3-2) 

Because of the finite width r ri , the mass, and thus the energy of the resonances, acquires an imaginary contribution. 
The loop diagrams of Fig. 1 then contain an imaginary piece, which describes the inelasticity of the NN scattering 
process above the pion threshold. The effect of direct pion production in NN scattering up to 1 GeV has been studied 
in Refs. [28, 47], and found to be very small in comparison to pion production via, e.g., the A-resonance. For the 
present purpose 7r (and r}) production will be restricted to occur only via resonance doorway states. 
We follow here the same parameterization of the penetrability as given in Ref. [33], i.e., 

T = jqv L {qR). (3.3) 

Here the momentum q is the cm. momentum of the irN system [or r)N system in case of the iV*(1535)], and should 
not be confused with the momenta in the NN system. The parameter R is fitted and chosen such that the product qR 
becomes dimensionless. We choose R = 6.3 GeV -1 for all three resonances under consideration. The dimensionless 
quantity 7 is the reduced width, and the penetrability factor vl is defined as 

VL{qR) = (qR)*nl(qR) + (qR)*jl(qRy M 

where L denotes the orbital angular momentum quantum number of the resonance. The functions jl and ul are the 
spherical Bessel and Neumann functions. Thus for the resonances A and TV* (1440) with the quantum numbers of P33 
and P\\ the width is given as 

q 3 R 2 

r A,fti = 7A,Pu 1 - ^ 2R2 &(q - qthr) (3.5) 

and for the resonance iV*(1535) with the quantum numbers Sn 

Tsu = ISuQ. 0(9 - qthr), (3.6) 

where the step function Q(q — qthr) ensures that the pion (eta) production starts at the corresponding thresholds 
in the NN cm. frame. The expressions for T ri show that the resonance width increases linearly with increasing 
momentum q and introduces a linear divergence in the imaginary part of the potential. Since this does not seem 
physical plausible, Ref. [33] introduced an additional phcnomenological form factor to the resonance width of the 
form 

2 g 2 K 2 J + K 2 

z{q } = mTr? ' (3J) 

so that Eq. (3.3) is modified to 

T = lq v L {qR)Z{q 2 ). (3.8) 



5 



As already mentioned in Ref. [33], the analytic form of the form factor is relatively unimportant, as long as the 
resonance curve is qualitatively described, i.e., it peaks at the resonance position. However, at the energies we are 
considering, setting Z(q 2 ) = 1 has only minor influence on the numerical calculation of the NN observablcs. Even so, 
we choose to employ Eq. (3.7) for the physical reason already mentioned. The form factor is defined such that it is 1 
at the exact resonance position, s = m^.., of the meson- nucleon system. The parameters 7 n are fitted such that the 
resonance width at the resonance position coincides with the experimentally extracted width [40] . A similar behavior 
has been observed in Ref. [30] where the resonance width is calculated as the imaginary part of the lowest order ANtt 
loop diagram. There a dipole cutoff with the same cutoff parameter A^A-n- as employed in the A A potential was 
used. 

Since we consider only single tt (or 77) production, we fit the width T ri to the corresponding partial widths as given 
in Ref. [40]. However, the A resonance decays with a probability > 99% into Ntt. Thus, in this case the partial width 
is equal to the total width, which lies between 115 and 125 MeV. For the A* (1440) resonance we use the partial 
width for the Ntt decay to determine the strength parameters. The total width has a wide ranging from 200 MeV to 
450 MeV Ref. [40]. We choose the value r%* = 391 MeV, from Ref. [48] to determine the partial width. 

The A* (1535) resonance decays mostly into Ntt and Nrj. Up to the ^-production threshold, which corresponds to 
the AA laboratory energy 1253 MeV, we only consider the partial width for Ntt decay, above the rj threshold we add 
also the partial width for the Nr] decay. The parameters we use for describing the width of the resonances included 
in our model are given in Table II. We tested the influence of the width cutoff parameter k on the AA observables. 
Varying k between about 200 and 1000 MeV did not lead to any visible difference in their description. Thus, for 
the A resonance we used the value (200 MeV) given in Ref. [33], while we for the two N* resonances we adopted an 
average value of 400 MeV. 

IV. FITTING PROCEDURE 

The parameters in our model are the masses of the mesons and baryons, together with the coupling constant 
and cutoff parameters for the baryon-baryon-meson vertices. With one exception the masses are fixed to the values 
recommended by the Particle Data Group (PDG) [40]. For the nucleon mass we employ the average of the neutron 
and proton mass, M = 938.926 MeV. Similarly, for the pion mass we use the average of the charged and neutral pion 
masses. The only mass parameters which are not fixed are those of the "er" mesons (one for each isospin channel), 
which we consider as free parameters. Since most mcson-baryon coupling constants are poorly known, we treat 
them as free parameters, which are allowed to vary within a certain range. We use the values of Refs. [30, 49] as 
starting points. The only exception is the NNtt coupling, which we fixed to <7^ rW7r /47r = 13.8, according to recent 
measurements at IUCF [50] and extractions by the Nijmegen group [51]. All cutoff parameters A are treated as free 
parameters and we allowed them to vary between 0.8 and 2.5 GeV. A more detailed discussion of the parameters 
obtained in the different fits will be given in the next section. 

Our goal is to find parameter sets for which our models give a good description of the A A observables, i.e., total 
and differential cross sections, as well as spin observables, up to 1.5 GeV. However, we do not fit directly to data but 
rather to a partial wave analysis, since in some cases the amount of data is rather sparse. As reference values we choose 
the SAID energy dependent phase shift solution Sp07 from the CNS Data Analysis Center [52]. However, we need 
to prepare from this solution a "data set" that is suitable for our fitting procedure. First, since the imaginary parts 
of our model are already fixed by fitting the widths of the resonances to their experimental values (see Section III), 
the partial wave inelasticities pi are already given. Therefore, we need to restrict ourselves to fitting only the phase 
shifts Si. Also, since we keep the pion coupling constant g%^^ fixed, we should not include higher partial waves in 
the fit. What constitutes a high partial wave is expected to depend on the energy considered. Therefore, we divide 
the energy range up to 1.5 GeV into three separate regions: up to 0.3 GeV laboratory kinetic energy we include only 
partial waves J < 2 as reference values, between 0.3 and 0.8 GeV we consider all partial waves with J < 4, and above 
0.8 GeV we include in addition all partial waves with J < 6. However, since there are only very few experimental data 
points for np scattering above 0.8 GeV, we exclude isospin-0 partial waves from the SAID data base above 0.8 GeV. 

In addition we need to account for the fact that the SAID energy-dependent solutions do not give errors. The 
absolute errors for each partial wave are estimated with a simple polynomial reflecting the errors of the single energy 
solution. In this way we try to account for the accuracy of the phase shift analysis, which varies with energy and 
partial waves due to different amount of available data. Especially partial waves with small absolute values get a very 
large relative error assigned, since they would be overrepresented if only simple relative errors were used. However, in 
order to avoid obscuring the plots, these errors are not included in our figures. In all cases, the assigned errors increase 
monotonically with energy and reach a maximum absolute value of 2 at 1.5 GeV. We tested different error setups 
and choose the one giving the best overall results. All partial waves get similar errors assigned. However additional 
weight-factors could be included during the fitting procedure itself, e.g. the 3 P2 partial wave above 800 MeV was 
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assigned a low weight in the final fits. After creating our "data base" according to these specifications, we perform a 
nonlinear least-square fit using the Minuit program provided by the CERN library [53] . A discussion of the resulting 
model parameters and the description of the phase shifts and observables will be given in the following section. 

The resonances are added one-by-one to our model, starting form a traditional meson-exchange NN scattering 
model. This way the influence of the different resonances can be understood separately from each other. 

V. RESULTS AND DISCUSSION 
A. Including the A(1232) - P 33 Resonance 

As a first step we consider only a meson exchange model together with the excitations of the A(1232) resonance. 
This is similar to the models used in Ref. [34, 49]. Here our focus is to create an updated model to which we can later 
add the Pn and Sn resonances. Thus, we first concentrate on the description of NN data up to 1 GeV, the energy 
regime which is expected to be dominated by the A(1232) resonance. In addition we include the contributions of the 
mesons ao, fo, and if, which all have a mass slightly below 1 GeV and were not included in [34, 49]. The influence 
these mesons have on the description of the phase shifts is relatively small and the coupling constants are determined 
to be of the order unity by our fitting procedure. They were allowed to vary within the interval [0.0,5.0] and the fit 
always stayed within these boundaries. 

The p-meson tensor-to- vector coupling ratio was allowed to vary within the limits of 5.9 and 6.3. The fit prefers a 
value of 5.9, to which the ratio was kept fixed in the subsequent investigations. In the case of the to meson this ratio 
is around -0.12 and of negligibly small influence. We thus fixed it to zero. We also studied the interplay between n 
and rf by alternately fixing one and letting the other one float. It turns out that, for the description of the NN phase 
shifts, only the combined strength of the two seem to matter, i.e., if one is fixed larger (on the order of 5) the other 
one will go to a smaller value, closer to 1. Therefore we always choose similar starting values for both of them, and 
let them vary between 1 and 6. Both couplings stayed within this range. Quark model prediction [54] give slightly 
higher values, but do predict both of them to be of the same magnitude. For the "<t" meson we let the coupling 
constants, masses, and cutoffs vary separately in the T = 1 and T = channels. However, the values for the masses 
always stayed relatively close together. 

Finally, we investigated the sensitivity of the A(1232) couplings to the overall description of the phase shifts. In 
[34, 49] the value of i/vATr/^ 7 *" = 0-35 is used. Similarly, the nN model from Ref. [55] employs f N ^/^ = 0.36. On 
the other hand, quark models relating f N ^ to f^NAir l ea d to values as low as 0.22 [56]. Because of this discrepancy 
we let the f N ^ n °at between those two values to see if the phase shifts show a preference. In addition we also let 
the cutoff A at att vary. The fit always prefers the lower value of the coupling constant and goes to the lower limit of 
the cutoff boundary, which was set for all couplings to 800 MeV. After lowering that bound even further, the cutoff 
Aatatt stabilizes at values between 600 and 700 MeV. Our fit runs consistently prefer a small value of Jnah- so we 
fixed this coupling for all further studies to 0.224. Since we independently fitted the width of the A resonance to 
its experimental value, this choice of coupling constant does not influence the size of the inelasticity provided by the 
A^A7r contribution. A similar investigation of the f N/ ^ p coupling did not reveal any sensitivity with respect to the NN 
phase shifts. Thus we kept the value for ffjAp/^ ^ ne same as in Refs. [34, 49] and only varied the cutoff parameter 

AwAp. 

In total we varied 22 parameters in this model, which includes the A(1232) resonance as the only means for pion 
production. It is well known that models of this kind are able to describe NN scattering up to about 1 GeV in 
semi-quantitative terms. A good parameter set obtained from our fit procedure is given as model (A) in Table III. 
Using those parameters we calculated pp and np observables as function of the projectile laboratory kinetic energy. 
The total cross section for elastic and inelastic pp scattering are shown in Fig. 3 as a long-dashed line together with 
the SAID Sp07 analysis [4] . As already pointed out in Ref. [34] , the inelastic contributions of a model of this type are 
not large enough to account for the experimentally observed ones. Consequently the total elastic cross section is being 
overprcdicted beyond 1 GeV. This tendency can also be observed in the differential cross section for pp scattering in 
Fig. 4, especially at the larger angles. A partial remedy would be double pion production, which starts to contribute 
above 590 MeV. We will not attempt to include such amplitudes in the present paper. The differential cross section 
for np scattering (Fig. 5) is experimentally less well determined than for pp and data above 800 MeV are sparse. 
Overprediction of the np differential cross section is only visible at the largest cm. angle. 

The analyzing power Am is shown in Fig. 6 for pp and in Fig. 7 for np scattering. Above about 800 MeV the 
analyzing power for pp scattering is strongly overprcdicted, especially for the larger angles. This phenomenon was 
already pointed out in Ref. [34]. The np analyzing power is also overprcdicted at similar angles, however at angles 
larger than 90° our calculation describes the existing data quite well. The spin correlation coefficients Ass, Ann, 
and All for different cm. angles are shown as function of projectile kinetic energy in Figs. 8-10 for pp scattering and 
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in Fig. 11 for np scattering. For the spin-correlation coefficient Ass the model predictions deviate from the data in a 
similar way as for the analyzing power, whereas the other spin-correlation coefficients are described reasonably well. 

In addition to the energy dependence of observables (for fixed angles) we also select three different fixed energies for 
cm. angular distributions of pp observables. Those energies are 400 MeV (Fig. 12), 800 MeV (Fig. 13), and 1300 MeV 
(Fig. 14) . As can be expected from the previous discussion, at 400 MeV our calculations are in reasonable agreement 
with the data, while deviations start to show up at 800 MeV and become more severe at 1300 MeV. 

The partial wave phase shifts of our model prediction (A) up to J=4 are shown as long-dashed lines in Fig. 15 
(T = 1) and Fig. 16 (T = 0). The imaginary parts of the partial wave amplitudes arc given in Fig. 17 (T = 1) and 
Fig. 18 (T = 0). 

B. Including the AT* (1440) - Pn Resonance 

The A* (1440) (Roper) resonance is the lowest J = \ nucleon excitation and thus the first nucleon resonance to add. 
Since it has the same quantum numbers as the nucleon this resonance couples to exactly the same meson fields as the 
nucleon, though not necessarily with the same strength. We are using the decay modes of the Roper as guide to which 
meson exchanges to include. It decays predominantly into the Nix channel (60-70%) [40]. A very crude estimate from 
the width of the Roper resonance indicates that the NN*ir coupling is about 40% of the NNtt coupling. In Ref. [57], 
Riska and Brown extract an average /^ rjv , 7r /47r = 0.0012 from the width and derive an even smaller quark model 
estimate. On the other hand, a coupled channel model for irN scattering [55] uses f^ N ,^/Air = 0.0024. We choose a 
value close to the latter, since the overall effect of the iterative AA*(1440) diagram is quite small. The second largest 
decay channel of the Roper is NAir (20-30%) and PDG gives a value of 5-10% for A(7T7r) s _ wavc . To mimic the latter 
effect , at least for the real part of the amplitudes, we introduce the "er" meson, and choose a relatively small value for 
its coupling. During the fitting procedure this value stayed small and none of our fits reached the large value favored 
in Ref. [55] . The calculations in Ref. [33] introduced additional diagrams containing p and to exchange in the iterative 
AW* (1440) diagram, though with relatively small coupling constants. As a test we included the oj meson, but since 
the overall effect in phase shifts and observables is minimal, we did not pursue the inclusion of vector mesons any 
further. 

Adding the A* (1440) resonance to model (A), without refitting the parameters, has only a moderate effect on both 
the T = 1 and T = phase shifts. However, to obtain a true assessment of the effect, the model parameters need to 
be readjusted. A qualitative fit to the phase shift including the A* (1440) resonance is given by the parameters listed 
in Table IV. The results are plotted as short dashed lines in Figs. 3-18. 

A major effect of the A* (1440) can be seen in the T — channels, which now exhibit pion production at threshold. 
Since AA contributions are forbidden, pion production close to threshold is provided solely by the decay of the Roper 
resonance. Overall, the total cross section for inelastic scattering is slightly increased. However, the contributions 
of the A* (1440) are by far too small to make up the difference with the experimental data. Though the changes in 
the x Pi and 3 P2 phase shifts are relatively large in going from model (A) to model (B), the final differences in the 
observables are much smaller. 

Thus, the inclusion of the Roper resonance has a small overall effect in the description of AA scattering. However, 
it is vital for the T = inelasticities close to threshold. Similar conclusions were drawn in Ref. [27]. We see in our fit 
that the pion-nucleon cutoff is slightly reduced when the Roper resonance is included. Thus we were tempted to see 
if we can force Annk to take a small value, e.g., 1 GeV, and still obtain a description of the data of equal quality by 
letting the NN*ir coupling take arbitrary values. In fact, it has been suggested [58] that a soft NN*ir form factor can 
be used in a boson-exchange model provided an additional pseudoscalar meson, with a mass slightly above 1 GeV, 
is introduced. We find, however, that despite the Roper resonance having the quantum numbers of the nucleon, the 
iterative AA* diagram with pion exchanges can not assume this particular role, and we can not get acceptable fits 
with a very soft pion form factor. 

C. Including the iV*(1535) — Sn Resonance 

The A* (1535) (Sn) resonance has the quantum numbers of the nucleon, except it is parity-odd. Thus the functional 
forms of the vertices have different spin structures, see Eqs. (2.12) and (2.13). The A*(1535) decays mainly into the 
Nir and Nr] channels, and thus we include only these mesons when considering the iterative diagram with the Sn 
resonance in the intermediate state. Values for the coupling g^ N ,^/Air vary considerably in the literature. The quark 
model estimate from Riska and Brown [57] suggests that it is slightly less then half of the NNtt coupling, while 
investigations of irN scattering give a value as low as 0.001 [59]. A study of r\ photoproduction from the deuteron 
[60] suggests a value of 0.1. Our value of 0.05 obtained from the fit is well within these boundaries. The AA data are 
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not particularly sensitive to this number as long as it stays around 0.1. For the coupling to the n meson we obtained 
0.99, which is higher than the value of 0.33 given in Refs. [59, 60]. 

When including the iV*(1535) resonance (without refitting), its effect is most prominent in the 1 Pi phase shift. In 
general, the T = phase shifts are more sensitive to the size of its contribution. However, a true assessment should 
only be made after a refit of the parameters. This gives the parameters of Table V, which constitute our model (C). 
The observables calculated with these parameters are plotted in Figs. 3-14, partial wave phase shifts in Figs. 15 and 
16, and the imaginary parts of the scattering amplitudes in Figs. 17 and 18 (solid lines in all figures). 

In general, the iV*(1535) resonance has a tendency to counterbalance the contributions of the iV*(1440) resonance. 
A description of NN observables with both N* resonances is of about the same quality as the description containing 
only the A resonance. Only a closer look at the inelasticities in the np phase shifts shows that these two resonances 
are necessary in order to obtain inelastic contributions in the region close to threshold. 

Comparing the parameters of model (C), given in Table V, with the parameters of the other two models shows that 
the coupling of the u) meson decreases when introducing the TV* (1535) resonance, indicating that this contribution 
has an overall repulsive effect. The p coupling, on the other hand, decreases minimally. In addition, the pion cutoff, 
A-nntt, assumes a slightly lower value. 

D. Role of the 3 P 2 Partial Wave in the Description of Polarization above 800 MeV 

Certain observables are particularly sensitive to specific phase shifts. For example, the analyzing power An is 
sensitive to spin-triplet states. Above roughly 800 MeV the analyzing power is strongly overpredicted for larger cm. 
angles, as seen in Fig. 6 for the pp excitation function. This effect is also seen in the corresponding np analyzing 
power, however not as strongly. Thus, one needs to concentrate on the large discrepancies seen in the T = 1 triplet 
phase shifts. The most obvious candidate is the 3 P2 phase shift, which above 800 MeV deviates strongly from the 
experimental analysis. The experimental analysis changes sign at about 1 GeV, while the model predictions remain 
positive and roughly constant. 

In order to isolate the influence of the 3 P2 partial wave on the observables, especially on the analyzing power, 
we replaced the phase-shift (5( 3 P2) from our model above 800 MeV with the one from the SAID analysis and then 
recalculated the observables. Since the inelasticity agrees with the SAID analysis, we stayed with our model prediction 
for this parameter. The resulting pp differential cross section and analyzing power are shown as function of laboratory 
kinetic energy for cm. angles 33° and 53° in Fig. 19, together with the original prediction for comparison. The large 
discrepancy almost completely vanishes. In Fig. 20 the differential cross section and analyzing power at 800 and 
1300 MeV are shown as function of cm. scattering angle. Again, once the 3 P2 partial wave is replaced with the SAID 
value, the discrepancy nearly disappears. All the other observables are only marginally affected by this change in 
the 3 P2 partial wave. The description of the analyzing power for np scattering at 53°, Fig. 7, is also improved by 
the change in the 3 P2 partial wave, while the backward angles are not affected and remain reasonably well described. 
For higher energies, the pp analyzing power is a smooth function of the energy at nearly all angles [4]. Thus we may 
speculate that once the fast change of sign of the 3 P2 partial wave between 800 and 1500 MeV wave is explained, the 
pp analyzing power can be described reasonably well. However, the present ingredients in our model can not account 
for a correct description of the 3 P2 partial wave. The efforts of Ref. [34] indicate that a variation of model parameters 
close to our model (A) is not able to improve the description of the analyzing power as function of energy. Instead, 
a different physical mechanism needs to come into play. A strong energy dependence in one single partial wave 
often indicates the opening of a new threshold in that particular channel, given by the coupling to another nucleon 
resonance. Since the 3 P2 partial wave has a very strong energy dependence above 800 MeV one may speculate that 
only a S'-wave coupling to this partial wave can be responsible, which would require a J = | resonance. The lowest 
J = | resonances have masses around 1.7 GeV and were found to be important for a description of irN scattering 
below 2 GeV [35, 61]. Consideration of either the Di 5 (1675) or the Fi 5 (1680) would be possible candidates to explore 
this avenue. 



VI. SUMMARY AND CONCLUSIONS 

In this paper we have studied the effects the first few resonances have on the description of the NN phase shifts 
in the energy regime around and above 1 GeV. We have incorporated the A(1232) resonance as well as the V*(1440) 
and iV*(1535) in a relativistic meson-exchange model. We started with a model including contributions of only the 
A(1232) resonance, which we refitted to NN data up to 1 GeV, and then added the contributions of the V*(1440) 
and V*(1535) resonances, now extending the fitted region up to 1.5 GeV. We verified the characteristic deficiencies in 
describing the pp analyzing powers above ^800 MeV that were pointed out already in Ref. [34] . We could isolate this 
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defect as being caused by a poor description of the 3 P2 partial wave above 800 MeV. This situation is not improved 
by including also the two N* resonances. We speculate that only a new channel, which couples with L = to the 3 P2 
partial wave will be able to cause such a strong energy dependence in that specific partial wave. Further investigations 
are necessary to investigate this assertion. 

We also find that, in general, the description of available np observables seems to be better than that of pp 
observables. However, this finding may be too optimistic due to the small amount of np observables above 1 GeV. 
We also observe that the total cross section for inelastic scattering is underpredicted by about a factor of two by 
our models, indicating that not enough inelastic channels are considered. In Ref. [34] this was compensated for by 
introducing a spin-independent optical potential providing an average absorption. As shown there, this provides 
a reasonable mechanism to incorporate a bulk contribution to account for single- and multiple-meson production. 
However, such an optical potential will not be able to account for the strong energy dependence in selected partial 
waves, e.g., the 3 P2 partial wave around 1 GeV. 

In this work the main effort was introducing contributions of two of the lowest N* resonances into the description 
of NN observables above 1 GeV laboratory kinetic energy. In order to have uniform description of the widths of the 
resonances we employ the relatively simple parameterization suggested in Ref. [33]. This can be improved by taking 
advantage of the recent work of the Excited Baryon Analysis Center (EBAC) [35] analyzing irN scattering in region 
below the 2 GeV nucleon resonance region, and e.g. employing vertex functions and self-energies determined in the 
■nN analysis in the NN system. This will allow to make connections between analyses of electromagnetic meson 
production reactions and the understanding of the NN system in this energy regime. 
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Threshold 


T lab [MeV] 


NNn 


285 


NA 


632 


AW* (1440) 


1136 


NNr] 


1256 


AA 


1355 


AW* (1535) 


1381 



TABLE I: Laboratory (kinetic) energy thresholds for single-pion and single-eta production in the AW system. In case of the 
resonances, the number indicates the laboratory kinetic energy for which the resonance becomes on-shell. 



Resonance 


Channel 


Tpartiai [MeV] 


R [GeV^ 1 ] 


7 


k [MeV] 


A 


Ntt 


~ 120 (a) 


6.3 


0.749 


200 


A"* (1440) 


Ntt 


~ 270 (t) 


6.3 


0.75 


400 


AT*(1535) 


Ntt 


~ 60 (a) 


6.3 


0.1286 


400 


AT* (1535) 


Nrj 


~ 79 (a) 


6.3 


0.4284 


400 



TABLE II: The parameters for the description of the partial resonance width as given in Eq. (3.8). The superscript (a) indicates 
values taken from Ref. [40], (b) indicates the value from Ref. [48]. 



Vertex 


a 


J p 


I 


m a 


'2 

9 NNa 
4tt 




A a [MeV] 


NNa 


TT 




o- 


1 


138.03 


13.8 




1519 




V 




0" 





547.3 


5.81 




830 




p 




1" 


1 


769. 


1.1 


5.9 


1281 








1" 





782.6 


23.1 


0. 


1382 




a(T = 


1) 


0+ 





497. 


5.505 




1807 




a(T = 


0) 


0+ 





480. 


3.545 




1990 




ao 




0+ 


1 


980. 


4.75 




1004 




rf 




0" 





958. 


1.62 




1433 




fo 




0+ 





980. 


2.5 




1274 










4-7T 


AjvAa [MeV] 


NAa 


IT 




0~ 


1 


138.03 


0.224 




640 




P 




1" 


1 


769. 


20.45 




1508 



TABLE III: Meson parameters of the model (A) including the A isobar as the only mechanism for pion production. J, P, and 
I denote spin, parity and isospin of the meson. Each meson vertex is multiplied by a form factor as given in Eq. (2.14), where 
n Q = 1 except for the NNp vertex where n a = 2. 



12 



Vertex 


a. 


J p 


I 


m a 


9 NNa 
4-7T 


f<x ha 


A a [MeV] 


NNa 


TV 




Q- 


1 


138.03 


13.8 




1498 




'1 




o~ 





547.3 


2.0 




1300 




n 
H 




i" 


1 


769. 


1.1 


5.9 


1281 




UJ 




1- 





782.6 


23.6 


0. 


1396 




a(T = 


1) 


0+ 





494. 


5.195 




1924 




a(T = 


0) 


0+ 





495. 


3.603 




1700 




do 




0+ 


1 


980. 


4.75 




1004 




rf 




0" 





958. 


1.24 




1433 




fo 




0+ 





980. 


2.5 




1471 










P 

NAq 
4-7T 


Ajvac [MeV] 


NAa 


TV 




0" 


1 


138.03 


0.224 




677 




P 




1" 


1 


769. 


20.45 




1508 










f' 2 , 


'i 


Ajvjv* q [MeV] 


AW*(1440)o< 


TV 




0" 


1 


138.03 


0.023 




800 




a 




0+ 


1 


494. 




0.74 


1119 



TABLE IV: Meson parameters of the model (B) including A isobar and the iV*(1440) as mechanism for pion production. J, P, 
and I denote spin, parity and isospin of the meson. Each meson vertex is multiplied with a form factor as given in Eq. (2.14), 
where n a — 1 except for the NNp vertex where n a — 2. 



Vertex 


a 


J p 


I 


m a 


'i 

9 NNa 


U/g a 


A Q [MeV] 


NNa 


TV 




Q- 


1 


138.03 


13.8 




1429 




V 




cr 





547.3 


2.0 




1298 




P 




i" 


1 


769. 


0.97 


5.9 


1310 








i" 





782.6 


20.5 


0. 


1421 




a(T = 


1) 


0+ 





503. 


5.254 




2000 




a(T = 


0) 


0+ 





530. 


4.164 




1516 




a 




0+ 


1 


980. 


4.75 




968 




v' 




<r 





958. 


1.57 




1424 




fo 




0+ 





980. 


2.31 




1442 










1 N Aa 
4tt 


AivAc [MeV] 


NAa 


TV 




0" 


1 


138.03 


0.224 




618 




P 




1" 


1 


769. 


20.45 




1422 










f' 1 . 

J JVJV*7T 

4tt 


'2 

9 NN*cr 


A NN * a [MeV] 


AW*(1440)a 


TV 




0" 


1 


138.03 


0.023 




800 




a 




0+ 





503. 




0.74 


1119 










y 

S NN*a 
4tt 


Ajvjv. q [MeV] 


AW*(1535)a 


TV 




0" 


1 


138.03 


0.05 




809 




V 




0" 





547.3. 


0.99 




962 



TABLE V: Meson parameters of the model (C) including A isobar and the nucleon resonances iV*(1440) and A r *(1535) as 
mechanism for pion production. J, P, and / denote spin, parity and isospin of the meson. Each meson vertex is multiplied 
with a form factor as given in Eq. (2.14), where n a — 1 except for the NNp vertex where n a = 2. 
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FIG. 1: The two-meson exchange iterative contributions involving nucleons (represented by solid lines) and nucleon resonances 
(double line). The dashed lines represent the exchanged mesons. The nucleon resonances and exchanged mesons are specified 
in the figure. 



1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I r ^ r 




q [MeV/c] 



FIG. 2: The width of the resonances A(1232) (solid line), iV*(1440) (dashed line) , and JV*(1535) (dash-dotted line) as function 
of the cm. momentum of the resonance. The dotted line indicates the width of the A resonance when no cutoff k is applied. 
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FIG. 3: The total cross section for elastic (left) and inelastic (right) pp scattering as function of the laboratory projectile kinetic 
energy. The solid line gives the prediction of the model C (Table V), containing the contributions of all three resonances P33, 
P11, and S11. The short-dashed line is based on the model B (Table IV) containing only the resonances P33 and Pn- The 
long-dashed line represents the results of model A (Table III) containing only the P33 resonance in addition to the one-meson 
contributions. The SAID analysis Sp07 [3] is indicated by the dotted line. The experimental data are referenced in the SAID 
database [4]. 




FIG. 4: The differential cm. cross section for pp scattering as function of the laboratory projectile kinetic energy at selected 
angles. The notation of the curves is the same as in Fig. 3. The experimental data are referenced in the SAID database [4]. 
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FIG. 5: The differential cross section for np scattering as function of the laboratory projectile kinetic energy at selected angles. 
The notation of the curves is the same as in Fig. 3. The experimental data are referenced in the SAID database [4]. 
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FIG. 6: The analyzing power An for pp scattering as function of the laboratory projectile kinetic energy at selected angles. 
The notation of the curves is the same as in Fig. 3. The experimental data are referenced in the SAID database [4]. 
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FIG. 7: The analyzing power Am for np scattering as function of the laboratory projectile kinetic energy at selected angles. 
The notation of the curves is the same as in Fig. 3. The experimental data are referenced in the SAID database [4]. 
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FIG. 8: The spin-correlation parameter Ass for pp scattering as function of the laboratory projectile kinetic energy at selected 
angles. The notation of the curves is the same as in Fig. 3. The experimental data are referenced in the SAID database [4]. 
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FIG. 9: The spin correlation parameter Ann for pp scattering as function of the laboratory projectile kinetic energy at selected 
angles. The notation of the curves is the same as in Fig. 3. The experimental data are referenced in the SAID database [4]. 
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FIG. 10: The spin correlation parameter All for pp scattering as function of the laboratory projectile kinetic energy at selected 
angles. The notation of the curves is the same as in Fig. 3. The experimental data are referenced in the SAID database [4]. 



21 




FIG. 11: The spin correlation parameters Ass, Ann, and All for np scattering as function of the laboratory projectile kinetic 
energy at selected angles. The notation of the curves is the same as in Fig. 3. The experimental data are referenced in the 
SAID database [4]. 
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FIG. 12: The angular distribution ol the differential cross section, the analyzing power and selected spin-correlation coefficients 
for pp scattering at 400 MeV laboratory kinetic energy. The notation of the curves is the same as in Fig. 3. The experimental 
data are referenced in the SAID database [4]. 
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FIG. 13: The angular distribution ol the differential cross section, the analyzing power and selected spin-correlation coefficients 
for pp scattering at 800 MeV laboratory kinetic energy. The notation of the curves is the same as in Fig. 3. The experimental 
data are referenced in the SAID database [4]. 
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FIG. 14: The angular distribution ol the differential cross section, the analyzing power and selected spin-correlation coefficients 
for pp scattering at 1300 MeV laboratory kinetic energy. The notation of the curves is the same as in Fig. 3. The experimental 
data are referenced in the SAID database [4]. 
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FIG. 15: The T = 1 NN partial wave phase shifts up to J= 4 as function of the projectile laboratory kinetic energy. The 
solid line gives the prediction of the model C (Table V), containing the contributions of the resonances P33, Pn, and Su. 
The short-dashed line is based on the model B (Table IV) containing only the resonances P33 and Pn. The long-dashed line 
represents the results of model A (Table III) containing only the P33 resonance in addition to the one-meson contributions. 
The SAID analysis Sp07 [3] is indicated by the dotted line. The dash-dotted line in selected phase shifts (here 3 F2 and 1 G4) 
represents the SAID analysis Sp04. 
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FIG. 16: The T = NN partial wave phase shifts up to J=4 as function of the projectile laboratory kinetic energy. The 
notation of the curves is the same as in Fig. 15. 
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FIG. 17: The imaginary parts of the T = 1 iViV partial wave amplitudes up to J=4 as function of the projectile laboratory 
kinetic energy. The notation of the curves is the same as in Fig. 15. 




FIG. 18: The imaginary parts of the T — NN partial wave amplitudes up to J=4 as function of the projectile laboratory 
kinetic energy. The notation of the curves is the same as in Fig. 15. 
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FIG. 19: The differential cross section (top panels) and the analyzing power An (bottom panels) for pp scattering as function 
of the laboratory projectile kinetic energy at two selected angles. The dashed lines are based on the calculation with the 
parameters of our model A (Table III). For the solid line the 3 P2 phase shift of the model is replaced by the SAID analysis. 
The experimental data are referenced in the SAID database [4]. 
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FIG. 20: The angular distribution of the differential cross section (left panels) and the analyzing power for pp scattering at 
800 MeV laboratory kinetic energy (top panels) and 1300 MeV (bottom panels). The notation of the curves is the same as in 
Fig. 19. The experimental data are referenced in the SAID database [4]. 



